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ABSTRACT 

We present new ATCA 21 -cm line observations of the neutral hydrogen in the nearby radio galaxy Centaurus A. We image in detail 
(with a resolution down to 7", ~ 100 pc) the distribution of Hi along the dust lane. Our data have better velocity resolution and better 
sensitivity than previous observations. The Hi extends for a total of ~ 15 kpc. The data, combined with a titled-ring model of the 
disk, allow to conclude that the kinematics of the H i is that of a regularly rotating, highly warped structure down to the nuclear scale. 
The parameters (in particular the inclination) of our model are somewhat different from some of the previously proposed models but 
consistent with what was recently derived from stellar light in a central ring. The model nicely describes also the morphology of the 
dust lane as observed with Spitzer There are no indications that large-scale anomalies in the kinematics exist that could be related 
to supplying material for the AGN. Large-scale radial motions do exist, but these are only present at larger radii (r > 6 kpc). This 
unsettled gas is mainly part of a tail/arm like structure. The relatively regular kinematics of the gas in this structure suggests that it 
is in the process of settling down into the main disk. The presence of this structure further supports the merger/interaction origin of 
the Hi in Cen A. From the structure and kinematics we estimate a timescale of 1.6 - 3.2 x 10** yr since the merging event. No bar 
structure is needed to describe the kinematics of the H i. The comparison of the timescale derived from the large-scale H i structure 
and those of the radio structure together with the relative regularity of the H i down to the sub-kpc regions does not suggest a one- 
to-one correspondence between the merger and the phase of radio activity. Interestingly, the radial motions of the outer regions are 
such that the projected velocities are redshifted compared to the regular orbits. This means that the blueshifted absorption discovered 
earlier and discussed in our previous paper cannot be caused by out-moving gas at large radius projected onto the centre. Therefore, 
the interpretation of the blueshifted absorption, together with at least a fraction of the redshifted nuclear absorption, as evidence for 
a regular inner disk, still holds. Finally, we also report the discovery of two unresolved clouds detected at 5.2 and 11 kpc away (in 
projection) from the H i disk. They are likely an other example of left-over of the merger that brought the H i gas. 

Key words, galaxies: active - galaxies: elliptical - galaxies: individual (Centaurus A) - galaxies: kinematics and dynamics - galaxies: 
structure - galaxies: ISM 



1. Introduction 

The assembly of early-type galaxies and its connection to the 
presence of an active nucleus (and in particular a radio-loud nu- 
, cleus) is still a matter of debate. Gas-rich galaxy mergers and 
i interactions could play an important role in providing the fuel 
' to trigger nuclear activity (including the radio phase). However, 
' recent studies of radio galaxies have shown that the activity in 
some of these galaxies may be associated instead with the slow 
accretion of (hot) gas. This is thought to be the case in particu- 
, lar, for edge-darkened radio galaxies where advective low effi- 
' ciency/rate flows could be the dominant mode of accretion (see 
e.g. Allen et al., 2006; Balmaverde et al., 2008). 

One way to address the issue of the relation between interac- 
tion/merger and AGN, is by studying the gas in radio galaxies. 
Nearby sources are prime targets as the morphology and kine- 
matics of the gas — both at large scales and in the nuclear re- 
gion — can be studied in detail and related to each other. Atomic 
hydrogen can be a useful tracer for these purposes. In particu- 
lar, H I gas found at kiloparsec scales, or larger, can be used to 
trace the evolution of the host galaxies (e.g. major merger vs. 
small accretions). On the other hand. Hi in absorption is one 
of the diagnostics to probe the central regions of active galactic 
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nuclei and the interaction of the AGN with its immediate en- 
vironment. H I absorption studies offer the unique possibility of 
traci ng the nature of the acc retion or the energetics of the outflow 
(e.g. lMorganti et al.Ll2005h . 

Unfortunately, the number of radio-loud galaxies where the 
atomic neutral hydrogen can be studied in such details both in 
emission and in absorption is limited. Due to its proximit}(3, 
the nearby radio galaxy Centaurus A (Cen A) is a unique ob- 
ject where to study the different phases of the gas, including Hi, 
in a radio-loud AGN. This can be done on spatial scales ranging 
from the inner kpc, where the influence of the AGN can still be 
relevant, up to many kpc where the signature of past merger(s) 
might be still visible through unsettled gas. 

Cen A harbours a well known r adio source classified as 
a Fanaroff-Riley type-I source (FR-I, lFanaroff&Ri"ievl [T97l 
with a total power Pi.vghz = 10^4 26 ^ Hz"'. The total ex- 
tent of the radio source reaches more than 8°, making Cen A 
the largest extragalactic radio source in the sky ("Coope r et al.L 
1965,) . The radio source shows a variety of features on dif- 
ferent scales. A radio jet extend s from sub-parsec scales (e.g. 
I Jones et al l lT99llHorhichi et al.L |2006() to a projected distance 



' At the assumed distance of Cen A of 3.8 Mpc dHarris et all 120091) , 
1 arcmin corresponds to 1.1 kpc. 
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of ~ 6 kpc from t he nucleus, ending in an inner radio lobe 
(iBums et al l ll983HClarke et all [l992^. Vast regions of fainter 
emission extend beyond the inner lobes: the middle lobe ex- 
tends out to 40 kpc, while the outer lobe is ev en more diffuse 
extending out to more than 500 kpc (see Fig. 1 in lMorganti et al.L 
[1991 . The different orientation and morphologies of the various 
structures could be the result of precession (e.g. iHavnes et al.L 
[l983) of the central engine combined with strong interaction 
with the surrounding medium. Alternatively, the nuclear activ- 
ity in Cen A has been recurrent, although recent spectral index 
studies suggest that the large-scale structure is currently under- 
going (or has very rece ntly undergone) some particle injection 
event (iHardcastle et al.L,20 09, and refs. therein). 

Most of the cold and warm gas in Ce n A is concen- 
trate d along the strongly warped dust lane (e.g. Nicholso n et al.L 
[1992: Ouillen et al... .2006.) . The neuti-al hydro gen present in 
this regio n (both in emission and a bsorption) has been stud - 
ied bv e.g.lvan der Hulst et al.l(ll983l) : lvan Gorkom et alJ(ll990l) : 
ISarma et al (12002 ). The morphology and kinematics of the outer 
regions of the H i disk suggest tha t some of the gas has no t yet 
settled into regular rotating orbits (Ivan Gorkom et al 1 [19901) . 

H I has also been found well outside thi s disk. At a distance 
betwee n 10 and 15 kpc from the nucleus, ISchiminovich et al] 
([T99I have found a partial ring structure with a smooth N-S 
velocity gradient and rotating perpendicular to the gas located 
along the dust lane (but co-rotating with the stellar body). In the 
north-eastern region of this ring, young stars have been detected 
that have possibly formed by the intera ction between the radio 
jet an d the neutral hydrogen in the ring dOosterloo & Morpantil 

Bool . 

Cen A was shaped through hierarchical merging over many 
Gyr as r evealed by the sprea d in ages of globular clusters 
(see e.g. IWoodlev et al.L l2010t). The last majo r merger event 
happened about 5 Gyr ago (" Peng et al.L 1200 4). The promi- 
nent warped dust and gas disk in the central region, and the 
outer partial ring structure, suggest that Cen A experienced a 
recent merging event where its l arge elliptical body merged 
with a smal le r gas-rich ga l axv dBaade & MinkowskiL Il954t 
iTubbs', '1980'). lOuillen et all dl993h calculated a timescale of 
about 200 million years since the core of the infalling galaxy 
merged with the centre of the elliptical galaxy. This is consis- 
tent with the timescales derived by the presence of tidal de- 
bris and by the shell-like features containi ng atomic hydrogen 
(Schiminovich et al., 1994:' Peng et al.Ll2002l) and molecular gas 
( iCharma nd aris et al., 2000). 

H I absorption is detecte d against t he nucleus, the northern jet 
and the southern lobe (e.g. 'Roberts, 1970; van der Hulst et al., 
ll983Hvan Gorkom et al.L[][990: Sai-ma et al. 2002). The absorp- 
tion against the jet and southern lobe has been interpreted as gas 
that is part of the large-scale disk w hich is ro ughly perpendic- 
ular to the jet axis (PA - 50°, Tin gav et all 1998). Since no 
absorption is found against the northern radio lobe, this lobe is 
likely pointing towards us whereas the southern radi o lobe p oints 
away. Van Gorkom et al. (1990) and iBurns et alJ (1 19831) con- 
cluded that the northern part of the gas disk is (at least partly) 
behind the radio source and the southern part of the H i disk is in 
front. For the component of H i absorption against the nucleus, 
earUer observations have shown the presence of redshif ted ab- 
sorption against the (unresolved) core (e.g . van der Hul st et all 
ll983Hvan Gorkom etal.'. 1990; Sarma et al., 2002), but our re- 
cent results, discussed in Morganti et al. (2008, Paper I) and in 
this paper, have shown that blueshifted absorption against the 
core is also present. 



Despite the many studies devoted to Cen A, there are a num- 
ber of open questions that are essential for the understanding of 
the formation and evolution of Cen A and its AGN. For exam- 
ple, the relation between the age of the radio structures and the 
timescale of the merger is not completely clear. In particular, it 
is unclear what triggers the current AGN activity (i.e. the activ- 
ity producing the inner radio lobe) and whether it is related to 
the merging event. Thus, it is important to obtain the morphol- 
ogy and kinematics of the gas disk to the largest possible radii 
and investigate the structure of the disk. This can provide an es- 
timate of the timescale of the merger and allow to investigate the 
presence and kinematics of unsettled gas (e.g. whether signifi- 
cant radial motions are present and, if so, on which scale). At 
the same time, it is important to investigate the distribution of 
H I down to the most inner parts of Cen A to understand whether 
kinematical signatures of gas connected with the fuelling of the 
AGN are present. 

The large amount of H i present in Cen A and its proximity 
allow to obtain deeper observations which can be used to inves- 
tigate the questions above. To achieve this, we have performed 
new spectral-line observations of Cen A that have higher spa- 
tial and velocity resolution and better sensitivity than those of 
previous studies. These new observations have already allowed 
us to detect blueshifted absorption against the nucleus while 
also the redshifted absorption extends to higher velocities com- 
pared to previously known. These results have been presented 
in Paper I. This gas has been interpreted as part of a circumnu- 
clear disk, likely the counterpart of what was already detected in 
CO and other lines, leaving still open the question of what trig- 
gers the current AGN activity. To further investigate this issue, 
we present in this paper the study of the kinematical state of the 
H I disk and the results from 3D modelling. 

The paper is organised in the following way. In Sect. 2 we 
present the new observations while in Sect. 3 we describe the 
results. In Sect. 4 we derive and discuss tilted-ring models of 
the Hi disk that are based on the 3D data cube that describe 
the strongly warped Hi disk as well as previously published 
data. We show that the Hi is dominated by rotation, with ra- 
dial motions giving a significant contribution in the outer parts. 
Furthermore, we argue that the blue- and redshifted absorption 
against the nucleus must be mainly caused by gas close to the 
nucleus. In Sect. 5 we discuss the consequences of the obser- 
vational and modelling results for the formation and evolution 
process of Cen A, the merger timescales involved and the im- 
plications for the AGN activity. A short summary is given in 
Sect. 6. 

2. Observations and data reduction 

Cen A was observed with the ATCA in April 2005. In Tab. [T] 
some basic information about the observations is summarised. 
These data were also presented in Morganti et al. (2008). Three 
different standard configurations (750 m, 1.5 km and 6 km) 
were used, each for a full 12-h run. This setup gives sufficient 
Mv-coverage as well as good spatial resolution. In addition, it 
avoids baselines shorter than 77 m. This is important because the 
very strong, very extended continuum emission of Cen A makes 
bandpass calibration in practice impossible on short baselines 
where the flux is well over 100 Jy (see below). 

The total bandwidth covered was 16 MHz - corresponding 
to a velocity range of 3000 km s"' - with 512 frequency chan- 
nels. This is much broader than used in ear lier observ ations (e.g . 
van der Hulst et all 1 19831: Ivan Gorkom etal , 1990; S arma et all 
|2002) . In order to have enough line-free channels for the contin- 
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Total on-source integration time 




3 X 12 h 


Observation dates 




tz/t4/ZZ Apill ZUUj 


Bandwidth 




16 MHz 


Number ol channels m the cubes 




512 


Channel spacing (dV) 




6.6 km s 


Velocity resolution (Av, after Hanning) 




1 J.z km s 


Shortest/longest baselme 




77 m/6 km 


Bandpass/flux calibrator 




DT/'C iC\1A iilO 


Gain/phase calibrator 
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Data cube 


high resolution 


low resolution 




LllliiUilll 


mniict (l 


Beam (HPBW) 


8.1" X 6.8" 


19.6" X 18.2" 


Beam position angle 
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-48.3° 


rms noise (mjy beam"') in the cubes 


1.3 


1.2 


rms noise (10" atoms cm"^) 


30.92 


4.45 


rms noise (Mq pc 
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0.36 


Peak flux continuum (Jy beam"') 
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Fig. 1 : H I emission in greyscale overlaid with the low resolution radio continuum (thin, grey contour levels: 0.5 Jy beam increasing 
with a factor of 1.5) and absorption contour levels on top (thick contours, red in the online version). Absorption contour levels: 
0.33 X lO'^ cm"^ increasing with a factor 3. Emission contour level: 1.5 X 10^° cm"^. 



uum subtraction, the central velocity of the observing band was 
about 200 km s"' below that of the systemic velocity of Cen A 
to account for the Galactic H i emission/absorption. The spec- 
tral resolution is 13.2 km s"' after Hanning smoothing. We used 
PKS 1934-638 as bandpass and flux calibrator. To monitor the 



system gain and phase changes, PKS 1315^6 was observed as 
a secondary calibrator for ten minutes every hour, although self 
calibration was also employed (see below). The data re duction 
was performed using the MIRIAD package (ISault et al.L[T995i) . 
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Fig. 2: Top panel: Total intensity contours over-plotted on DSS image. Contour levels: 0.2, 0.4, 0.8, 1.6 (absorption only), 3.2, 6.4, 
12.8, 25.6 and 51.2 x 10^'^ atoms cm"^ for the emission (black) and absorption (white). Tspin = 100 K. The horizontal line indicates 
the position of the pv-slices in Fig. |3]and[8] Bottom left panel: spectrum of H i absorption cloud. Bottom middle panel: spectrum of 
the nuclear absorption. Bottom right panel: zoom-in of bottom middle panel. 



The extremely strong radio continuum requires care- 
ful bandpass calibration in order to obtain flat spectral 
baselines. Following a tec hnique applied successfully by 
lOosterloo & Morgantil ( 12005 ^ on earlier observations of Cen A, 
we smoothed the bandpass calibration obtained from PKS 1934- 
638 with a box-car filter 15 channels wide. Since the instrumen- 
tal features in the bandpass are fairly broad in frequency, this 
procedure is allowed and it minimises the contribution of extra- 
noise from the bandpass calibration that otherwise would have 
occurred since the bandpass calibrator is fainter than the very 
extended emission of Cen A. Through this smoothing, the effec- 
tive flux level of PKS 1934-638 increases to about 60 Jy which 
is higher than the detected flux on almost all baselines. The 
smoothed bandpass correction was applied to the unsmoothed 
data of Cen A. This procedure ensures that the velocity resolu- 
tion is preserved without increasing too much the noise of the 
data. The strongest H i absorption is just over 2 Jy and the dy- 
namic range in the spectral line cubes is better than 1:1000. 

The continuum was subtracted using UVLIN by making a 
second-order fit to the line-free channels of each visibility record 
and subtract this fit from the spectrum. The line cube, as well as 
the obtained continuum image, was cleaned. Banning smooth- 
ing was applied to the line cube. Frequency independent self- 
calibration of the phase using the channel with the strongest nu- 
clear absorption improved the image quality, but a few channels 
with the strongest absorption still show some calibration arti- 



facts. The final cube was derived by combining the datasets of 
the three observations. The beam size in the resulting uniformly 
weighted cube is 8'.'1 x 6'.'8, PA = -12°3. The image noise is 
about 1.3 mJy beam"'. For the kinematic modelling of the disk 
(see Sect. 5), a cube with robust weighting set to was made 
while also tapering the visibility data with a Gaussian taper. The 
resulting beam size for this cube is 19'.'6 x 18'.'2, PA = -48°3. 
This lower resolution cube is a compromise between spatial res- 
olution and sensitivity. 



3. Results 

3. 1 . The structure of the H i emission 

As expected from previous observations, H i is detected both in 
emission and in absorption and most of the gas is concentrated 
in a relatively edge-on, warped structure that mainly follows the 
dust lane (see Fig.[T]and|2]i. On both the eastern and western side, 
the H I extends beyond the visible dust lane and, in total, has an 
extent of ~ 15 kpc in diameter. The total mass of the H i disk is 
4.9 X 10^ Mq. This is slightly higher than what was derived from 
VLA observations (4.2 x 10^ M q, after correcting for th e dif- 
ference in the assumed distance; Ivan Gorkom et allll990 ). The 
difference is likely due to better separation of emission and ab- 
sorption in the central regions (due to our better spatial resolu- 
tion), and to the higher sensitivity of our observations. 
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In Sect. 4 we present detailed modelling of the geometry and 
the kinematics of the H i, but we briefly discuss the main points 
here. The main conclusion is that the kinematics of the H i is reg- 
ular in the inner regions of the disk. Large-scale radial motions 
exist, but are only present at larger radii. There are no indications 
of anomalies in the kinematics of the gas (on the scale down to 
the inner beam, ~ 100 pc) that could be related to supplying 
material for the AGN. 

Interestingly, the radial motions of the outer regions are such 
that the projected velocities are redshifted compared to the regu- 
lar orbits (see below and Fig. [3). This means that the blueshifted 
absorption discovered in Paper I (see also Sect. 3.2) cannot be 
caused by out-moving gas at large radius projected onto the cen- 
tre. Therefore, the interpretation of the blueshifted absorption, 
together with the already known redshifted nuclear absorption, 
as evidence for a regular inner disk, can still hold. 

Previous studies have already shown that the kinematics of 
the gas in the dust lane is dominated by rotation (see ;?v-diagram 
in Fig. 2 of Paper I). Although there is a fair degree of symme- 
try in the H i distribution, significant deviations from symmetry 
are present, not only in morphology, but also in kinematics. As 
can be seen from Fig.|2] the famous optical dust lane has a clear 
counterpart in Hi. Below the dust lane, south-west of the nu- 
cleus, a similar amount of H i i s present as in the dust lan e it- 
self. Previous observations (e.g. ISchiminovich et all [r994 l) had 
already shown the presence of this Hi, but it is better visible 
in our data. Most of the H i south-west of the centre can be ex- 
plained in a simple way. In the optical, the dust lane is only vis- 
ible where the gas and dust disk is in front of the galaxy, i.e. the 
north-eastern side. South-west of the nucleus, the disk is behind 
most of the stellar body and it will be much less visible in ab- 
sorption. However, the H i of the disk is seen in emission and 
therefore both the back and the front side of the disk are visible 
equally well. 

The kinematics of the H i also shows that the outer gas disk 
south of the centre is not perfectly regular. The position-velocity 
plot given in Fig. [3] shows the kinematics of part of the south- 
western disk. Since the position-velocity plot was taken far from 
the nucleus, both the emission and absorption should reflect the 
kinematics of the same outer H i disk. The only difference be- 
tween absorption and emission should be that the absorption 
comes from that part of the south-western disk that happens to 
be in front of the radio continuum lobe, but the kinematics of the 
absorption should smoothly fit that of the emission. However, 
the absorption is offset in velocity with respect to the emission, 
something that is not expected for a regular rotating disk with 
circular orbits. In addition, the absorption structure has a faint 
continuation in emission, see the box in Fig. [3] 

The presence of this gas is interesting for various reasons. 
The part in emission represents a spatially and kinematically 
continuous structure that appears to be, at least partly, connected 
to the large-scale disk. The gas originating the emission and the 
absorption (inside the box in Fig. [3]) is likely to be one contin- 
uous structure, representing front and back of one continuous 
structure (called "outer arm"). The observed broad profiles sug- 
gest that radial motions must be affecting the gas in addition to 
the rotation. Further intriguing is the fact that we do not detect 
similar structures on the opposite (northern) side of the dust lane. 
This already points to an asymmetry at large radius in the distri- 
bution of gas that is clearly a complication in the modelUng of 
the data (see Sect. 4). 



Table 2: Absorption-line properties. 



Component 


V 




T 


A'Hi/^s 




[Km/sj 


[Km/sj 








(1) 


(2) 


(3) 


(4) 


Nucleus 


550 


18 


0.83 


2.1 A X 10'" 


Nucleus 


574 


23 


0.62 


2.61 X 10'" 


Nucleus 


592 


20 


0.26 


9.52 X 10"* 


Nucleus (shallow) 




400 


* 


* 


Jet 


550 


19 


0.64 


2.23 X 10" 


Jet 


544 


25 


0.64 


2.93 X 10" 


Cloud 


454 


13 


0.07 


1.72 X 10"* 



(1) Velocity of the peak. (2) FWHM (3) Peak optical depth. (4) 
Hi column density [cm"- K"']- 
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Fig. 4: Velocity channel of the high-resolution cube (greyscale). 
The white contours show the Hi absorption (-15, -12, -9, -6 
and -3 cr). The thick black contours give the CO emission lLisztl 
(2001). In thin black contours (0.5, 1.0, 2.0 and 4.0 Jy beam"') 
the low-resolution continuum image is shown. The beam of the 
H I cube is indicated in the lower right corner. 



3.2. Hi in absorption 

The ATCA observations reveal new features also for the H i ab- 
sorption. Figure [T] shows the total intensity of the Hi in absorp- 
tion as well as the low-resolution continuum as obtained from 
the line-free channels. The radio continuum shows the inner 
lobe s and the inner jet, in agreement with p revious observations 
(e.g. Clarke et al., 1992-. Bu rns etali ll983). The Hi absorption 
is detected against the southern radio lobe, the northern jet and 
the unresolved nucleus, see Fig.|2]for more details. In addition, 
we detect absorption up to ~ 20" (corresponding to ~ 3 beams) 
east of the core, in the direction perpendicular to the jet (Fig.|4]l, 
indicating the presence of background radio continuum extended 
in that direction. There is a close correspond ence between this 
absorption and the CO emission as detected bv lLisztl(l200lh rsee 
Fig.©. 

Figure |5] shows a position-velocity plot along the jet PA. 
This illustrates how complex the absorption is and, in particu- 
lar, that there is redshifted absorption everywhere. The outer fil- 
ament at large radius is redshifted and partly in front of nucleus 
(see Sect. 3.1 and Fig.[3]l. It is therefore difficult to disentangle 
whether at least part of the redshifted absorption is at all associ- 
ated to infalling gas. No equivalent, large-scale blueshifted gas 
seems to be present. 
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Fig. 3: Position- velocity plot obtained along a slit in PA = 90° as illustrated in Fig.|2] The data is integrated over 24" in Declination 
to better highlight the faint outer H i arm structure in emission (black, inside the white box) which connects with the absorption 
(white) on the western side. Note that the greyscale is inverted (i.e. absorption white, emission black) compared to Fig.|4l|5]and|8] 
to better highlight the faint emission. 
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Fig. 5: Position-velocity plot obtained along the radio jet, as il- 
lustrated in the inset in the lower right corner. Contour levels: -3 
(white) and 3cr (black). The systemic velocity is indicated by the 
dashed line. 



As shown in Fig. |2] the absorption against the northern jet 
and the beginning southern lobe has a typical column density 
(assuming Tspi„ - 100 K) of the order of 10^' cm"^ which is 
of the same order as the column densities of the neighbourhing 
emission. Because both seem to be connected to the large-scale 
disk (Sect. 4) this suggests that 100 K is a reasonable choice for 
the spin temperature in this region. However, the south-western 
parts of the absorption have more than an order of magnitude 



lower column densities - i.e. < 10^° atoms cm"^ (see Fig.|2]l - 
and hence these column densities are below the detection limit 
(3cr level per channel) of the Hi emission. 



While the extended Hi absorption against the northern jet 
and the southern lobe is in agreement with previous observa- 
tions, most exciting is the result that the absorption against 
the unresolved nucleus is much broader than previously known 
(Fig.|2]i. In particular, the blueshifted part of the absorption was 
undetected in previous studies. A complete analysis and discus- 
sion of the origin and nature of this absorption against the nu- 
cleus has been presented in Paper I. While the shallow part of the 
nuclear absorption is new, the deep c omponents close to t he sys - 
temic velocity are in agreement with Ivan der Hulst et al ] (ll983h . 
Peak velocities, width of the profiles, optical depths and column 
densities are summarized in Tab.|2] 



For the nuclear components of the H i absorption, the value 
of the spin temperature is much more uncertain as the gas could 
be located closer to the AGN (see Sect. 4 and 5) and therefore 
T^ n\n = 100 K cou ld underestimate the real value. According 
to ' Bahcall&EkersI d 19691) and iMalonev et al] ( Il996h the spin 
temperature can easily reach a few times 10^ K for gas close 
to an AGN. Assuming a spin temperature of a few 1000 K, we 
estimate column densities of ~ 10^^ - 10^^ cm"^, similar to 



what is found in other radio galaxies (e.g. Peck & Tavloil 120011: 



iMomjian et an. l2002'.'2003':'Morganti et alil2004 <). We can com- 
pare the H I column density of Cen A with the value derived from 
X-ray observations. IE vans et al.l (120041) have estimated a column 
density for the nucleus (from Chandra and XMM-Newton data) 
that is comparable to our value if we as- 
3000 K. It is worth noting, that the column den- 
sity derived for Cen A from the X-ray is relatively high com- 
pared to other ra dio galaxies of similar radio power (see e.g. 
iBalmaverde etall 120081) . 



of 10^^ atoms cm ^ 
sume Tspin 
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Right Ascension I 

Fig. 6: H I emission obtained from the ATCA observations (greyscale; orang e in online version ) superimposed to an optical BVR 
image obtained after unsharp masking and adaptive histogram equalisation bv lPeng et alj (|2002|) . 



3.3. Two newly discovered H i clouds in the environment of 
CenA 

H I at large distances from the dust lane was detected and stud- 
ied by [Schiminovich et al. (1994). These large-scale structures 
are also detected in our data as illustrated in Fig. |6] which shows 
an overlay of the H i emission on an optical image obtained 
by iPeng et all (l2002h . The optical image emphasises the low- 
contrast features and in particular the complex and extended 
set of shells and the faint dust extensions. However, our obser- 
vations include only one pointing on the centre of Cen A and 
therefore these structures are not imaged as well as in the work 
of Schiminovich. We therefore consider the study of this large- 
scale H I structure beyond the scope of this paper. 

However, it is interesting to note that two additional unre- 
solved clouds are detected in our observations, far away from 
the Hi disk (see Fig. [1]). The first cloud is detected in absorp- 
tion against the NE radio-lobe (see also Fig. |2]i at a distance of 
4'.7 (about 5.2 kpc) from the nucleus (PA - 50°5) and has a ve- 
locity width (FWHM) of only Av = 19 km s ^ The velocity of 



this cloud (454 km s" ) is below the systemic velocity of Cen A. 
The optical depth is r = 0.05 corresponding to a column density 
of about 10^" atoms cm"^, resulting in an upper detection limit 
of 10^ Mq (assuming a spin temperature of 100 K and a filling 
factor of unity). The second cloud is located 1 1 kpc from the nu- 
cleus at a position angle of about 48°5 w.r.t. the nucleus and is 
detected in emission at v = 418 km s"'. The cloud is only visi- 
ble in the low resolution cube, which indicates a diffuse nature. 
Also this cloud has a small velocity width (Av = 24 km s"') 
and an Hi mass of about 8 x 10^' Mq. Both clouds are located 
along the jet axis where also ionised gas has been dete cted (e.g. 
iDufour & van den Berghi[l978HMorganti et al.Lll99lh . 

4. The H I model 

The modelling of the Hi has at least two main goals. It allows to 
investigate the formation and evolution of Cen A (e.g. its merger 
history) through the study of the morphology and kinematics of 
the gas over the whole disk. The molecular and ionised gas disks 
are not as extended as the H i. The second reason is to pin down 
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Fig. 7: Parameters of the H i model described in Sect. 4. 



the location of the gas responsible for the nuclear absorption. 
Is the absorbing gas located in the vicinity of the nucleus or at 
large distance projected in front of the core? The answer to this 
question is important for the understanding of the AGN activity 
(see Sect. 1 and Paper I). 

A number of modelling studies of the large-scale gas disk of 
Cen A have been done in the past (see Morganti , 2010^ for a re- 
view ), but previous H i emission-line studies ([van Gorkom et al.l 
ll99Ql) did not allow such a detailed investigation as presented 



here. lBland et alj(ll987h and lNicholson et all d 19921) used the ve- 
locity field derived from Ha observations to obtain models de- 
scribing the warped dust-lane structure with a set of tilted-rings. 
The largest radius fitted corresponds to about 245" (~ 4.5 kpc 
see also Fig.fTOll. Ouillen et al.l (11992) obtained from CO obser- 
vations a first warped disk model that matches the optical appear- 
ance of the dust lane and which is also cons i sten t with the H i ro- 
tation curve derived by van Gorkom et al.l (1199 0). In the most 
recent model of Ouillen et al. (2006) derived from new mid- 
infrared observations, the position angle is quite similar — at 
large radii — to the one of i Nicholson et al ] (ll992h although the 
inclination is closer to edge-on for r > 140" (see also Fig.fTOli. 

ISparkd d 19961) developed a number of dynamical models to 
explain the geometry of the disk as that of a near-polar struc- 
ture precessing around the symmetry axis of an approximately 
oblate galaxy potential. T hese models compare w ell with the ori- 
entation of the model by INicholson et all ( 11992 !). According to 
all these models, the gas disk is strongly warped and crosses 
the line-of-sight several times at different radii. With this com- 
plex geometry, it is impossible to unambiguously assign veloc- 
ities to a position on the sky and, therefore, to derive a reli- 
able tilted-ring model ("Rogs tad & Shost^ 11972 ') and a rotation 
curve based on the velocity field only. Thus, a construction of 
model cubes is needed for a proper comparison with the data. 

To achieve this we have created full 3D model cubes us- 
ing GASMOD, a phyton code to calculate the tilted-ring models 
(kindly provided by Paolo Serra). The initial parameters of in- 
clination and posit ion angle for the large-scale disk were taken 
from Ou illen et al.l {2006). Because these parameters do not ex- 
tend to the entire range of radii covered by the H i, we have ex- 
trapolated them to larger radii. A flat rotation curve, that rises 
steeply inside our first ring has been assumed. The best values 
for the coordinates of the centre and for the systemic velocity re- 
sulted in agreement with the value from literature values. For the 
systemic velocity the value of Vsys = 546 km s ' was used, sim- 
ilar to the value of 'van Gorkom e t aP (11990 . v'sy, = 542kms-i). 
All parameters (surface brightness, rotation velocity, inclination, 
position angle) were then independently varied from their initial 
values for every single ring. 
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Fig. 9: Left panel: Spitzer image (lOuillen etalil2006h over plotted with intensity contours of the high resolution H i model. Right 
panel: Total intensity map of the H i model at the resolution of the Spitzer data. 



The model cubes were compared with the Hi data cube to 
make sure that the model would describe the kinematics of the 
observed gas. The changes to the model were made based on a 
visual inspection of model cube and data. We also produced from 
the model the distribution of the total intensity of the Hi that 
was then compared with the Spitzer image in order to make sure 
that the model would also reproduce the overall characteristics 
of this image. In order to perform this comparison, models were 
produced with the spatial resolution of the Spitzer data (see be- 
low). 

The starting parameters clearly did not provide a satisfac- 
tory description of the data. The final value of the parame- 
ters as function of the radius are shown in Fig. |7] The rota- 
tion curve is essentially flat and midly decreasing with radius, 
a p henomenon also observed in other early-type galaxies (see 
e.g. iNoordermeer et all l2007h . Figure |8] shows a comparison 
between a slice of the data and the starting and final model. 
The slice is close to what is shown in Fig. [3] and, therefore, it 
includes some of the unsettled gas described in Sect. 3.1 that 
we do not expect to see described by the model. However, the 
region of strong Hi emission is clearly better reproduced by 
our Hi model (figu re on the left) compared to the model of 
lOuillen et al.1 ( l2006l) . In addition, our Hi model also describes 
the dust morphology (Spitzer image) as shown in Fig.|9] 

Thus, the model (Fig. |7]i describes most of the Hi mor- 
phology and kinematics confirming the regular rotation of the 
Hi structure for r < 6 kpc. The orientation parameters (position 
angle and inclination) confirm the heavily warped structure of 
the disk. In addition, the model shows that in H i a gap in the dis- 
tribution of the gas is present between 10" and ~ 45" (from 
0.2 to 0.8 k pc), confirming \ yhat yyas already found by other 
studies (e.g. Nicholso n et al.L Il992t lOuillen et all 1 19921 l2006t 
[Espada et al., 2009) for the dust, molecular and ionised gas (e.g. 
Ha and CO). 

Our best model uses parameters different from what was 
found in previous modelling. Figure [TO] shows a comparison be- 
tween the par ameters of our best model and the parameters from 
lOuillen et al.l (12006i) . The main noticeable difference is in the in- 



clination that appears to stay closer to edge-on for r > 90" in our 
model. Remarkable is the similarity of our parameters with the 
independent results from the study of the stellar light emerging 
from the kiloparsec-scale, ring-like structure detected using the 
VLT dKainulainen et al. J2009l) . 

Our model description of the data is solely based on cir- 
cular rotation and does not require a bar st ructure (previously 
claimed to be present bv lMirabel et al.L 'l999). in agreement with 
the results from the Spitzer images of Keene_et al. (2004) and 
lOuillen et all (I2006h. 

Recently, 'Es pada et al.l (l2009l) have modelled high- 
resolution CO data obtained with the SMA for the central kpc of 
Cen A and invoked the presence of a weak bar in these central 
regions. Interestingly, while their model seems to very poorly 
describe the jump in velocity observed in the CO velocity field 
at about 15" from the centre (see their Fig. 4 compared with 
their Fig. 11), our model naturally explains this as only due to 
the geometry of the disk. The presence of a s mall-scale bar is 
also ru led out by the near-IR data presented bv lNeumaver et all 
(120071) . 

5. Discussion 

The new data presented here have allowed to image in detail the 
H I emission and absorption of the warped disk of Cen A and to 
identify new features in the gas. In this section we discuss the 
implications of the H i results for the formation and evolution 
(e.g. merger history) of Cen A and the different phases of AGN 
activity. 

5.1. Modelling implications: The large-scale H i distribution 

Our data confirm that the structure and kinematics of the H i in 
Cen A is dominated by rotation, down to the central regions, and 
most of the H i (but clearly not all) is confined to a fairly regular, 
warped disk (for r < 6 kpc). 

In addition to this, we detect outer unsettled gas that appears 
to be located in an arm-like structure. Therefore, some of the 
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gas in emission and absorption is observed below the dust-lane 
and it is not reproduced by our modelling. This arm shows a 
smooth velocity gradient - similar for absorption and emission 
but with an offset in velocity between the two. Because of the 
regular velocity characteristics of this structure, it is likely that it 
will become part of the outer part of the rotating disk in the near 
future. 

This type of structure is expected as a result of a merger 
or accretion as shown in many numerical simulations (e.g. 
[Hibbard & van Gorkom, 1996; Barnes, 2002). Thus, the Hi fur- 
ther support the hypothesis of Cen A being the result of a rela- 
tively recent merger of an elliptica l galaxy with a smaller gas- 
rich (SMC-like?) d isk galaxy (e.g. iBaade & Minkowski! Il954t 
lOuillenet all 11993 . see also Sect. 1). 

Using the results of our modelling, we can derive a rough 
estimate of how long ago the merger must have taken place to 
allow a large fraction of the Hi to settle down in the warped 
disk we observe. We can assume that the gas located at larger 
radii than we can model (i.e. r > 350") has performed at least 
one or two revolutions since the merger took place giving a 
rough estimate resulting in 1.6 - 3.2 x 10^ yr. T his time esti- 
mate i s in agreement, inside the uncertainties, with Ouillen et al.l 
(1 1 9931 2 x 10^ yr) and Peng et al. (2002, 3 x 10** yr), while it 
is shorter than what derived by Sparke ( 1996) - 7.5 x 10^ yr - 
that would correspond to almost five revolutions. Such a long 



timescale would be sufficient to relax most of the gas in the disk, 
in contradiction with our H i results. 

The merger hypothesis does fin d further support fr o m the 
partial ring structure discovered by ISchiminovich et all ( 11994 
see also Sect. 1). In that respect the existence of additional, 
low mass (< 10'' M©) Hi clouds is not surprising. The ve- 
locities of the two newly discovered Hi clouds (see Sect. 3.3) 
are in rough agreement with the spatial velocity distribution 
of the "outer ring" structure. Clouds and filaments of ionised 
gas have been found in t his region al ong the jet direc tion (e.g. 
Duf our & van den Berghl [1978: Mor ganti etdl Il99lh . For the 
origin of these structures, two hypothesis have been put for- 
ward. In addition to the merger origin, an alternative scenario 
was the cooling condensations in a tenuous, thermally unsta- 
ble hot component of the interstellar medium gas. This scenario 
was somewhat favoured to explain the filaments of ionised gas, 
as it would explain th e highly turbulent velo city structure ob- 
served there (see also lOosterloo & Morgantii 12005) . However, 
such high velocity gradients are not seen in the H i profile of the 
newly detected clouds which could mean that the clouds were 
not hit by the jet. Whether the newly discovered Hi clouds are 
a cooling condensation or are more likely the left over of the 
interaction/merger discussed above, remains unclear. 

As a final remark, we would like to note that also the distri- 
bution of the H I confirms what was found by other studies of the 
molecular gas and dust (Nicholson et al., 1992; Ouillen "et al.L 
2006), i.e. the lack of gas in region between 10" and 45" (0.2 kpc 
to 0.8 kpc) while on the smallest scales circumnuclear disks are 
found (Paper 1 and ref. therein). The peak in the surface bright- 
ness around 150" from the centre suggests a ring-like structure 
in the distribution of H i in Cen A (see Fig.|7]). H i ring-like struc- 
tures are often detected around early-type galaxies and there are 
a number of possible scenarios for the formation of such a struc- 
ture (see e.g. Donovan et al. 2009 for a summary and references). 
A ring morphology can be the result of the presence of a bar 
structure but we found no evidence for such a structure in Cen A 
and indeed many other galaxies with gas rings do not show bars. 
In some cases, "burning out" of the gas in the central regions 
of the ga laxies has been p roposed, but the famous example of 
IC 2006 dSchweizer et all [T989[) shows that also gas accretion 
can form these structures. 



5.2. Merger and different phases of radio activity 

The kinematics of the H i in Cen A is dominated by rotation 
down to the nuclear regions. The discovery that the nuclear ab- 
sorption is also partly blueshifted, (see Paper I) has questioned 
the interpretation that H i against the nucleus represent evidence 
for the direct fuelling of the AGN. Furthermore, our analysis 
presented here shows that the nuclear blueshifted absorption is 
not due to gas at large radius seen in projection against the core. 
The nuclear absorption is likely originating from a nuclear disk 
that represents the H i counterpart of the structure seen from the 
molecular gas (see Paper I for a detailed discussion). This nu- 
clear disk is part of the overall warped disk structure. Thus, our 
results confirm what suggested in Paper I that the nuclear H i gas 
does not constitute direct evidence of gas infall into the AGN. 

However, it is worth mentioning that radial motions (in ad- 
dition to rotation) on the scale of tens of pc ha ve been detected 
in the high resolution, near-IR data obtained bv lNeumaver et all 
(l2007l) . Surprisingly they found that for higher excitation lines 
[SiVI] and [Fell] the velocity pattern is increasingly dominated 
by a non-rotating component, elongated along the radio jet. 
Interestingly, these non-rotational motions were detected along 
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the jet with redshifted velocities (compared to the systemic) 
seen on the main-jet side and blueshifted on the counter-jet side. 
These motions (stronger in [Si VI]) can be explained as backflow 
of gas accelerated by the plasma jet. This gas can perhaps be in- 
volved in the fuelling of the AGN. Whether a similar situation is 
also occurring in the Hi gas can only be investigated with deep 
VLBI observations. 

In order to investigate in more detail whether a link can be 
found between the merger/accretion event and the nuclear ac- 
tivity, it is important to compare the timescales of these events. 
Cen A shows a complex radio structure with lobes on different 
scales, from a few kpc of the inner lobe (see Fig. [T]) to tens of 
kpc for the middle lobe and up to several hundred kpc for the 
outer lobe (Sect. 1). Th ese three radio struct ures could be either 
the result of precession (iHavnes et all 1 19831 see also Sect. 1) or 
the nuclear activity in Cen A has been recurrent. However, the 
large-scale structure is currently experiencing a particle i njection 
event, as recent spectral-index studies suggest ( Hardcastle et al.L 
l2009l and refs. therein). 

The a ges of these lobes are not easy to estimate. According 
to ISaxton et al. (2001) the age of the northern middle lobe can 
be estimated to be about 1 .4 x 10*^ yr assuming this lobe has been 
created by an old episode of jet activity that has been interrupted 
by the disruption of the jet by the merger. Since then, this struc- 
ture has been rising buoyantly away from the nucleus. If this is 
the case, the merger that has brought the H i would actually be 
the cause of the disruption of the old radio activity. The timescale 
of the merger derived from the H i could support this idea. The 
time scale of the inner radio lobe, ~ 5-10x1 O^yr dCroston et al.L 
12009 . and refs. therein), would indicate that it takes quite some 
time for the system to recover from this disrupting event and 
for the gas to settle again in a steady accretion flow. A (possi- 
ble) time-delay between the merger and the onset of the AGN 
has also been detected in a few other nearby radio galaxies (e.g . 
Tadhunter et al., 1 9961: lEmonts et all 120061: iDavies et all 120071: 
Wild etal.. .20101 

It is worth n oting that if the scenario proposed by 
ISaxton et al.l (12001 ) is correct, the delay between the disruptive 
merger and the onset of a new phase of activity could be related 
to the time re quire d by the X-ray halo to build up again (see 
ISansom et all 1200 0). 

On the other hand, recent spectral-index studies have sug- 
gested different timescales for the large-scale radio structure. 
High frequency studies show a spectral break for the south- 
ern lobe that would s uggest a life time of 3 x 10^ yrs (see 
iHardcastle et"al]. |2009|) . Furthermore, no spectral break is ob- 
served in the northern lobe (both the giant and the middle lobe), 
suggesting that the y have v ery recently undergone some parti- 
cle injection event d Alvarez et al.. 2000; H ardcastle et al., 2009, 
and refs. therein). Thus, if the timescales derived from the spec- 
tral index are correct, it would indicate that that the overall radio 
emission starts relatively late compared to the merger. 

In both cases, the comparison of the timescale derived from 
the H I and the radio structure, together with the relatively reg- 
ularity of the H I down to the sub-kpc regions, do not suggest a 
one-to-one correspondence between the merger and the phase of 
radio activity. In particular, if the timescales derived from spec- 
tral index are correct, more than one radio burst can be produced 
while the merger is settling in. This would indicate that, while 
the merger could be responsible for bringing gas in the vicinity 
of the nucleus, the accretion mechanism that fuels the AGN is 
likely not related to the merger. 

This is in agreement with the idea that in radio galaxies of 
relatively low radio power (i.e. FR-I type) the fuelling of the 



black hole proceeds directly from the hot phase of the interstel- 
lar medium in a manner analogous to the Bondi accretion. As 
sug gested by a number of authors (see e.g. Allen et al., 200^ 
iHardcastle et all 120071) the hot. X-ray emitting phase of the in- 
tergalactic medium (IGM) is sufficient to power the jets of sev- 
eral nearby, low-power radio galaxies. In the case of Cen A, a 
mass accretion rate M Bondi = 6.4 x lO^'^M© yr"' and a Bondi 
efficienc y of ~ 0-2 % has been derived from Chandra obser- 
vations dEvans et all l2004h . Howev er, it would decreases i f the 
new valu e of the BH ma ss derived bv' Neumaver et all d2007h and 
Cappella ri et al.l d2009h is used. These values are consistent with 
what is found for other FR-I radio galaxies dBalmaverde et all 
2008) while the efficiency appears to be lower than the "canoni- 
cal" 0. 1 for the high efficiency optically thick accretion disks. On 
the other hand, and to complicate the picture, it is known (from 
the Fe Ka line and the l arge column densi ty measured in the X- 
ray towards the nucleus !e vans et aTll2004l) that Cen A has large 
quantities of cold gas in a molecular torus at about Ipc from the 
black hole, more typical of FR-II galaxies. 



6. Summary 

In this paper we have reported ATCA H i synthesis observations 
which have higher spatial and velocity resolution than previous 
studies of atomic hydrogen in Cen A. Our study concentrates 
mainly on the distribution and kinematics of the gas associated 
with the dust lane of this galaxy. These data show that the kine- 
matics of the H I in Cen A is dominated by rotation down to the 
nuclear regions. 

The warped morphology and kinematics of the Hi can be 
described with a tilted-ring model with gas on circular orbits. 
Regular kinematics of the gas is observed for r < 6 kpc. We 
find no need for non-circular motions between 1-6 kpc. There is 
a lack of Hi between 0.1-1 kpc, confirming results from ob- 
servations of CO, ionised gas and dust. The parameters of the 
model (in particular the i nclination) are different from what de- 
rived in previous models dOuillen et all 120061 : iNicholson et all 
1992) but nicely c onsistent with those from the study of the cen- 
tral stellar ring of iKainulainen et al] d200^ . Outside 6 kpc ra- 
dius, asymmetries and arm structures are observed as well as 
radial motions. 

The Hi fraction in Cen A (Mjj, /Lb - 0.01) is relatively 
low for an early-type galaxy, indicating that the H i disk might 
well be the result of accretion of a small (SMC-like?) galaxy. 
Based on the regular kinematics we suggest that this accretion 
occured about 1 .6 - 3.2 x 10^ yr ago, in agreement with previous 
estimates. This is too old to directly trigger the current phase of 
AGN activity (the age is a few x 10^ yr for the inner radio lobes), 
but might have interrupted the previous episode of AGN activity. 

We find extensive H i absorption, against nucleus, jet and 
lobes. Part of the absorption (in particular with velocities close to 
Vsys) is part of the regular rotating disk, while some of the absorp- 
tion against the southern lobe is part of an outer arm structure. 
A large fraction of the nuclear absorption is not produced by 
outer gas confirming our hypothesis from Paper I. This is in par- 
ticular the case for the newly discovered blueshifted absorption. 
This absorption is, at least partly, due to a circumnuclear disk at 
r < 100 pc. Whether part of the (redshifted) gas is infalling and 
fuelling the AGN is difficult to establish with our data. The re- 
sults of iNeumaver et al.] (1200 7^ — who detected radial motions 
along the jet axis in near-IR observations — suggest that it could 
be worth exploring with deep VLBI data whether such a compo- 
nent also exists in H i. 
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All in all these results suggest that in galaxies with an AGN 
where a merger/accretion event has clearly occurred in the re- 
cent past, the connection between these two events (merger and 
activity) is far from clear. In the case of Cen A, the AGN may in- 
stead be fuelled via a low efficiency/rate accretion, as suggested 
for other FR 1 galaxies of similar radio power. The merger may 
actually have played the opposite role of temporarily disrupting 
the radio jet and producing the northern middle radio lobe. The 
delay in creating the much younger inner radio lobe, may be con- 
nected to the time that it takes the gas to settle again in a steady 
accretion flow. 
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